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Abstract. We have computed optical Zeeman spectra of magnetic white dwarfs for field strengths between 10 and 
200 MG and effective temperatures between 8000 and 40 000 K. They form a database containing 20 628 sets of 
flux and circular polarization spectra. A least-squares optimization code based on an evolutionary strategy can 
recover relatively complex magnetic field topologies from phase-resolved synthetic Zeeman spectra of rotating 
magnetic white dwarfs. We consider dipole and quadrupole components which are non-aligned and shifted off- 
centre. The model geometries include stars with a single high-field spot and with two spots separated by ~90°. 
The accuracy of the recovered field structure increases with the signal-to-noise ratio of the input spectra and 
is significantly improved if circular polarization spectra are included in addition to flux spectra. We discuss the 
strategies proposed so far to unravel the field geometries of magnetic white dwarfs. 
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1. Introduction 



About I 3% of all white dwarfs ha ve strong magnetic 



fields 



method devised for the analysis of r apidly rotati n g mag- 
netically active main sequenc e stars ( Semel 198S ; Donati 



jeLween 10 s and 10 q Gauss flWickramashighc fe 



Ferrario 2000| ; [Jordan 200 In many of these magnetic 
white dwarfs (MWDs), the surface field geometries devi- 
ate from simple centred dipoles. This holds for isolated 
MWDs and for the MWDs in accreting close binaries 
( Wickramasinghc fc Ferrario 2000; Bchwope 1995). While 



higher modes are often thought to decay on timescales 
longer than the t <; 10 10 yr of the fundamental mode, 
Muslimov et aL (1995) showed that quadrupole or oc- 
tupolc components may survive via the Hall effect if 



ct al. 1989| ; |Brown et al. 199l|) is not applicable to MWDs. 
The field geometry of MWDs can be derived, however, 
from the analysis of the pure Zeeman splitting of the pho- 
tospheric lines and their circular polarization properties 
as a function of rotational phase. Because of the large 
Zeeman shifts, this approach must include the whole op- 
tical range for B <; 50 MG. In the absence of positional 
information from the Doppler effect, however, the inver- 
sion of the flux and polarization spectra is an intricate 
task. Trial-and-error fits of centred or shifted dipoles and 



an intqnial Luiuidal magnetic field cumpuiiuiiL is pieseul. 
Therefore, studies uf Lhu surface field sLiucLme piuvide 

clues on the internal field configuration and its influence 
on the evolution of MWDs. 

The photospheric spectra of hydrogen-rich MWDs are 
characterized by broad absorption structures formed by 
the overlap of numerous components of the Balmer lines, 
shifted by hundreds or even thousands of A from their 
zero-field positions by the linear and quadratic Zeeman 
effects. These shifts completely obliterate the Doppler 
shifts caused by rotation even in the most rapidly rotating 
MWDs. As a consequence, the Zeeman-Doppler imaging 



quadrupoles ( Wickramasinghe & Cropper 1985; Putney & 



Jordan 1995) are incapable of exploring the full parameter 



space of possible solutions. We present a new strategy us- 
ing a pre-computed database of synthetic MWD spectra 
and an automatic quality-of-fit optimization algorithm. 

A first approach along these lines was presented by 
Donati et aL (1994), who used a maximum entropy algo- 
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rithm (MEM) to fit a matrix of areal filling factors for a 
grid of synthetic flux and circular polarization spectra to 
simulated input data. This way, they obtained the 'sim- 
plest' and, according to Occam's razor, most likely fre- 
quency distribution of transverse and longitudinal field 
strengths over the visible stellar disc at each rotational 
phase, a so-called ZEBRA plot, but this approach does 
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not provide any information about the spatial structure expand into four equations (Beckers 1969; Hardorp et al 



of the nagnetic field. 



In this paper, we investigate to what extent the un- 



1976) which can be solved by one of several different al 
gorithms: (a) the method of Wickramasinghe & Martin 



der lying global magnetic field distribution can be recov- 
ered directly from least-squares fits to phase-resolved flux 
and polarization spectra of a given signal-to-noise ratio. 
This approach uses the spatial information provided by 
the magnetic fields seen at different rotational phases and 
has the advantage that the uncertainties in the parame- 
ters describing the global field structure can be directly 
related to the noise in the spectra. Its disadvantage lies 
in the necessary restriction to model fields which can be 
described by a sufficiently small number of parameters. 

We assume fields which can be represented by cen- 
tred or offset dipole and quadrupole components which 
need not be aligned with each other. The specific geome- 
tries tested here include a star with a single high-field 
spot and one with two spots separated by ~90°. Our com- 
puter code allows us to calculate areal filling factor ma- 



(1979D assumes that the source function is linear in the op- 



trices analogous to ZEBRA plots, the resulting flux spec- 
tra, and the wavelength-dependent circular polarization 
for a given magnetic field model viewed at a number of 
rotational phases. We compare the results with the ref- 
erence input (which are simulated data in this case) and 
determine the best-fit parameters using an evolutionary 
optimization strategy. 



tical depth and that between two successive depth points 
the Stokes parameters can be described by exponential 
functions; (b) direct Runge-Kutta integration; (c) acceler- 
ated A iterations (Takeda 1991); (d) an approximation for 
large Faraday rotation (Ramaty 1969); or (e) matrix ex- 
ponential solutions ( pittmann 1995| ) Intensive tests per- 
formed by H. Schmidt and S. Jordan in Kiel have demon- 
strated the numerical equivalence of these methods with 
high accuracy. For the present paper, we have calculated 
an extensive database of synthetic flux and circular polar- 
ization spectra using Ramaty's approximation, which is 
always justified in white dwarf atmospheres, and is rather 
efficient with regard to CPU time. 

The temperature and pressure structure of our atmo- 
spheres is t aken from zero field LTE models ( [Koestei 



et al. 1979 ). As a conseq uence, the m agnetic pressure 
and magnetic blanketing ( lordan 1992| ) have been ne- 



glected. Convection is assumed to be suppressed by the 
field (lordan 2001). For the line opacities, data from the 
Tubingen group (Forster et al. 1984; Rosner et al. 1984; 
Wunncr et al. 1985) were used. Bound- free opacities were 
calculated using a modified approximation by Lamb & 



tra computed tor a wide range ot ncld strengths, viewing 
angles, and effective temperatures. Section ^ describes the 
general design of the magnetic field models and Sect. ^ the 
construction of the integrated spectra from the database 
for a given model of the magnetic field. Section explains 
the optimization code, describes the specific field models 
subjected to the reconstruction tests, and investigates the 
ability of the code to deduce the respective field param- 
eters from the integrated flux and polarization spectra. 
Finally, the power and also the limitations of our approach 
are discussed in Sect. ||. 

In forthcoming papers, we will analyse phase-resolved 
spectral flux and circular polarization data of MWDs ob- 
tained at the ESQ VLT. 



2. The database 

2.1. Radiative transfer for Magnetic White Dwarfs 

Our synthetic Zeeman spectra and wavelength-dependent 
polarization data are computed with the most recent ver- 
sion of the code developed by S. Jordan. The polarization 
originates from the different absorption coefficients n\, K T , 
and k p for left- and right-handed circularly polarized light, 
and linearly polarized light travelling perpendicularly to 
the magnetic field, respectively, and is described by the 
four Stokes parameters /, Q, V, and U. The influence of 
the Faraday rotation and the Voigt effect is accounted 
for by the magneto-optical parameters pr and pw- The 



Thf present paper is arranged as toiiows. m sect, y Sutherland (1974Q which leads to small errors only Jordan 



we describe tne database ot nux and polarization spec- 1992 | ; [Merani et al.T995|; |Jordan fc Merani 1995|) "and 



saves an enormous amount of computing time. The ap- 
proach described here was developed in two diploma the- 
ses ( Euchne^W9^_ Rjdml_999 ) and was also implemented 
by [Burleigh et al. (1999j ). 



three radiative transport equations of Unno (1956) then 



2.2. Database spectra 

We computed a three-dimensional grid of Stokes I and 

V model spectra with the effective atmospheric temper- 
ature T, the magnetic field strength B, and the field di- 
rection ijj relative to the line of sight as the independent 
variables. We considered T = 8000 K, 9000 K, 10 000 K, 
11 000 K, 12 000 K, 13 000 K, 15 000 K, 17000 K, 20 000 K, 
25 000 K, 30 000 K, 40 000 K, B = 10 MG to 200 MG in 
steps of 1 MG, and ip = 0°, 29°, 41°, 51°, 60°, 68°, 76°, 
82°, 90°, i.e., equidistant in cos?/'- This yields a database 
containing 12 x 191 x 9 = 20 628 model spectra for J and 

V each. All spectra are calculated for a surface gravity of 
logg = 8. Since we do not include the linear polarization, 
the field direction is sufficiently constrained by the total 
field strength and the longitudinal field component. As an 
example, Fig. [j] shows a section around Ha for a sample of 
database spectra with T = 15 000 K, B = 20 MG and five 
angles of ip, equally spaced in cos"0- A typical property 
of these Zeeman spectra is the weak angular dependence 
of Stokes I, except near 0° and 90°, and the more pro- 
nounced dependence of Stokes V. Somewhat simplified, 
Stokes I carries much of the information on the distribu- 
tion of the absolute value of B over the surface of the star, 
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Fig. 1. Examples of the database spectra for B = 20 MG, T = 15 000 K, and five angles of ip equally spaced in cos-0 
covering the Ha cr~-, 7T-, and er + -components. a) Intensity, b) degree of circular polarization. The spectra are shifted 
vertically by arbitrary amounts to avoid overlap. 



while Stokes V is needed to derive the distribution of the 
field directions. 

If limb darkening is important the direction cosine \i 
of the line of sight with respect to the vertical direction in 
the stellar atmosphere needs to be considered as a further 
parameter in the database. Hence, including a wavelength- 
dependent description of limb darkening requires an ex- 
pansion of the number of model spectra in the database 
by a factor equal to the number of /z-values considered. 
For the present calculations, we use a simple limb dark- 
ening law which is independent of wavelength and avoid 
this extension of the database (see Sect. 4.3 below). 



3. Magnetic field geometry 

A curl-free field which originates only from sources in the 
stellar interior can be described by a multipole expan- 
sion of the scalar magnetic potential, using spherical har- 
monics with coefficients I and m = 0, . . . , I, which describe 
the zonal and sectoral periodicity of the field, respectively 



(see, e.g., Langel 1987). The number of free parameters 
of the field geometry is 1(1 + 2), i.e. 8 (15) for I = 2 (3). 
Defining the viewing geometry requires three additional 
parameters for the orientation of the rotation axis and 
the inclination. 

The optimization procedure adopted for the present 
tests can handle only a limited number of multipole pa- 
rameters and becomes inefficient already when the oc- 
tupolcs (I = 3) are included. We have restricted the com- 
plexity of the field, therefore, by including only the two 
lowest zonal harmonics, commonly referred to as 'dipole' 
(1 = 1, m = 0) and 'quadrupole' (I = 2, m = 0), allowing 



their axes to be inclined with respect to each other. We 
do not consider the m = 1 and m — 2 quadrupoles, but 
instead include a common offset of the dipole-quadrupole 
combination from the centre of the star. This hybrid model 
has ten free parameters: two polar field strengths; two an- 
gles each for the directions of the axes relative to the ro- 
tation axis; the three components of the offset; and the 
inclination, i.e. the angle between the rotation axis and 
the line of sight. An offset from the centre was included 
because of its popularity and simplicity (e.g., the Earth's 
magnetic field is approximately that of a shifted dipole). 
The chosen field structure deliberately includes some very 
similar field geometries described by different sets of pa- 
rameters: a combination of aligned dipole and quadrupole 
can be approximated by a shift of the dipole. At sufficient 
signal to noise, the reconstruction procedure can distin- 
guish between such geometries, a result which is of inter- 
est by itself. While our hybrid model is useful for the tests 
performed in this paper, its limited complexity may not 
suffice for the interpretation of real, observed spectra. 

We consider rotating MWDs viewed at an inclination 
i with respect to the rotation axis. Note that a frac- 
tion / = 0.5 (1 — sin i) of the stellar surface is perma- 
nently hidden from view and that phase-resolved Zeeman 
spectroscopy provides no information on the field on this 
hidden fraction of the surface. In order to save comput- 
ing time, we restrict ourselves to simultaneously fitting 
flux and polarization spectra at four rotational phases, 
<j) = 0.0, 0.25, 0.5, and 0.75. We avoid a special geometry 
by choosing = not to coincide with the nearest ap- 
proach of one of the axes to the line of sight. 
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Fig. 2. a) Rotational geometry of the MWD models. The x-axis points towards the observer, the zo-axis marks the 
rotation axis. The dashed line marks the intersection of the x-z-plane with the a;o-yo _ pl ane - i denotes the inclination, 
and u> the rotational phase angle, b) Magnetic geometry of the MWD models. The z'- and z"-axes mark the axes of 
symmetry for the dipole and quadrupole components. The lines of intersection of the x'-z' -plane and the x"-z"-plane 
with the xo-yo-plane are given by the dashed and the dotted lines, respectively. 



Since observational restrictions often prevent taking 
phase-resolved data, we also consider the amount of in- 
formation which can be retrieved from a single flux and 
polarization spectrum or even a single flux spectrum only. 
In this case, the data provide information on the magnetic 
field structure only for one hemisphere of the star. 

At any given phase 0, the polarization depends on the 
components of the field transverse and parallel to the line 
of sight. In order to describe these components, we intro- 
duce four Cartesian coordinate systems (Figs. ||a and||b): 
(i, ii) systems £' and in which z' and z" describe the 
dipolar and quadrupolar axes of symmetry, respectively; 
(iii) the observer's system S, in which the x-axis points to- 
wards the observer and the z-axis lies in the plane defined 
by the x-axis and the rotation axis; and (iv) the auxiliary 
system So with zq the direction of the rotation axis which 
defines the inclination angle i. The rotational phase an- 
gle ui = 2n<t> is defined with respect to the direction of the 
Xo-axis which lies in the x-z-plane for lj = 0. 

The components of the surface field B(r') of the cen- 
tred dipole in system £' are 

{B%, =3S p d ol xV/(2r' 5 ), (1) 

(B%=W* ol y>z>/(2r> 5 ), (2) 

(B d ) z ,=B« ol (3z> 2 ~r> 2 )/(2r> 5 ), (3) 

with r 1 = [x' , y', z') and \r'\ 2 — r' 2 = x' 2 + y' 2 + z' 2 . Cor- 
respondingly, the components of the centred quadrupole 
in S" are 



(B q )*« 


= B« ol x"(5z" 2 


-r" 2 )/(2r" 7 ), 


(4) 


(£<V 


-B; ol y"(5z" 2 


-r" 2 )/(2r" 7 ), 


(5) 


(B q ) z " 


= B^z"(5z" 2 


-3r" 2 )/(2/' 7 ), 


(6) 



with r" = (x", y", z") and |r"| 2 = r" 2 = x" 2 + y" 2 + z" 2 . 
£' and S" are tilted with respect to the rotation axis zq 
by angles d and 9 q , respectively. The azimuth angles 
of the tilt in So are $ d and $ q at phase = 0. We apply 
the appropriate rotations and the translation to transform 
the fields into the observer's system E, add the dipole 
and quadrupole components, and obtain B(r) for each 
surface element, with B x = B\ = Bcostp the longitudinal 
component of the field. 

The angle 5 between the dipole and quadrupole axes 
and the angles r/ d and i] q between the line of sight at 
phase <p and the dipole and quadrupole axis, respectively, 
are given by 

cos 5 = cos e d cos 8 q + sin 9 d sin 6 q cos($ q - $ d ) , (7) 
cos r] d = cos i cos 9 d + sin i sin d cos(27r^ + <& d ) , (8) 
cos rf 1 = cos i cos q + sin i sin q cos(27r</> + $ q ) . (9) 

Our magnetic geometries were selected for the purpose 
of providing sufficient complexity for an effective test of 
our reconstruction routine. The offset r' g from the centre 
is a simple means of producing a substantial amount of az- 
imuthal asymmetry if r' D s is perpendicular to the dipole 
axis, while r' Q g parallel to the dipole axis allows to test the 
ability of the routine to distinguish between aligned cen- 
tred dipole-quadrupole combinations and a shifted dipole. 

4. Input spectra for the reconstruction procedure 

4.1. Integration of the database models 

We divide the spherical star into surface elements defined 
by equal steps in latitude and longitude. For given distri- 
butions of the magnetic field vector B and the effective 
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Fig. 3. Spectral flux (top) and circular polarization (bot- 
tom) for centred dipoles viewed pole-on with polar field 
strengths of (A) 200 MG, (B) 80 MG, and (C) 30 MG for 
T = 15 000 K. The library-based spectra computed using a 
database of four /x- values (grey) are compared to those for 
\i = 1, corrected with a mean limb darkening law (black) . 
Spectra (A) and (B) have been shifted upwards to avoid 
overlap (1.1 flux units each). The insert shows the relative 
flux differences for the 30 MG case (C). 

temperature T over the surface, let a be the running in- 
dex of the surface elements which are visible at a given 
rotational phase > (f> > 1 and which have sizes A a , cen- 
tral field strengths B a , field directions ip a , and direction 
cosines fx a . The Stokes parameters and (Vx) are then 
computed as weighted sums of the individual contribu- 
tions corrected for limb darkening by a factor (dis- 
cussed below) 



(h) 
(Vx) 



(<t>) 



E 

a(4>) 



LD 



(10) 



We represent the wavelength-dependent contributions I\^ a 
and Vx >a from surface element a by appropriate interpola- 
tion in the database grids of the parameters T, B, and ip. 
For T and ip, a bilinear interpolation suffices. For the field 
strength, we consider all spectra representative of the B- 
variation over the finite surface element. We found that a 
number of 900 surface elements per hemisphere is a good 
compromise between CPU time and needed accuracy. This 



number is sufficient to avoid spectral structure caused by 
the finite clement size. 

Our code can account for temperature variations over 
the surface of the white dwarf, but in this paper we con- 
sider only stars with uniform surface temperatures. 

4.2. Relation to the observed flux 

The flux observed from a star of radius R at distance d is 



fx = Fx R 2 /d 2 



(11) 



where F\ = n (I\) is a function of T, for a given magnetic 
field distribution. The interpretation of observed Zeeman 
spectra in terms of (I\) and (Vx) involves T and R/d as 
fit parameters. 

For the present tests, T and R/d are considered as 
fixed parameters and the quantities fitted by variation of 
the field parameters are Fx and V/I — (Vx) / 

4.3. Limb darkening 

We have compared (i) the full radiative transfer for each 
surface element which accounts for the /i-dependence and 
the variation of B across the element already in the atomic 
data, (ii) an interpolation between the spectra for dis- 
cretized // and B. and (iii) the application of a wavelength- 
independent linear limb darkening law replacing the in- 
terpolation in fi. Method (ii) uses spectra calculated for 
/i = 0.1, 0.4, 0.7, and 1.0. Method (iii) employs a linear 
law with coefficients which are valid for the visible wave- 
length range, 



= I x (fi)/h^=i = 0.70 + 0.30 M ■ 



(12) 



Test calculations for the three approaches (i) to (iii) were 
performed for centred dipoles viewed pole-on with polar 
field strengths of (A) 200 MG, (B) 80 MG, and (C) 30 MG. 
We found the differences between (i) and (ii) to be minute. 
Case (iii) differs by a wavelength-dependent factor which 
arises from the neglect of any wavelength dependence in 
the limb darkening approximation. Figure ^| shows the re- 
sults for cases (ii) and (iii) at an effective temperature of 
15 000 K. For all three field strengths, the spectra com- 
puted for case (iii) deviate by at most 5% from those of 
cases (i) and (ii). The insert shows that the absorption 
lines are about 2% deeper than for the correct treatment. 

The simple limb darkening law of Eq. [l^ is en- 
tirely acceptable for the present tests which interpret 
synthetic spectra with spectra of the same origin. The 
above comparison suggests, moreover, that a wavelength- 
independent linear limb darkening law may also be ac- 
ceptable for the interpretation of observed optical Zeeman 
spectra of white dwarfs. 

5. Reconstruction of the field geometry 

In this Section, we describe a variety of magnetic field 
and viewing geometries and test the ability of our code 
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Fig. 4. Model (Al), centred dipole viewed at i = 60°: a) B-tp diagram, b) field strength and longitudinal component, 
c) flux and polarization spectra. The + and — symbols indicate the sign of the longitudinal component of the magnetic 
field. For clarity, the flux spectra at <j> = 0.25, 0.5, and 0.75 have been shifted upwards by one flux unit each. 



to reconstruct their parameters from flux and circular 
polarization spectra at <f> = 0.0, 0.25, 0.5, and 0.75. All 
calculations were performed for an effective temperature 
of 15 000 K. In order to simulate real data, noise was 
added to the input spectra at the four phases as de- 
scribed below. Because of the added noise, the recon- 
structed field is not necessarily identical to the input field. 
For the present tests, the wavelength range was restricted 
to 4000 < A < 7600 A, which contains the most important 
Balmer line components, and all spectra were rebinned 
into 10 A bins, yielding 361 data pixels per spectrum, a 
total of 1444 pixels in the combined flux spectra at four 
rotational phases, and another 1444 pixels in the polariza- 
tion spectra. 

5.1. Magnetic field models 

We define seven different magnetic field and viewing ge- 
ometries against which we test our reconstruction code. 
The field configurations (A) to (F) are characterized by an 
increasing level of complexity. The geometrical and spec- 
tral properties of the models are summarized in Figs. ^-0- 
In each case, the centre panel shows the distributions of 
the total field strength B and of the absolute value of the 
longitudinal component B\ over the visible hemisphere at 
the four selected phases. The + and — symbols indicate 
the sign of the longitudinal component. The range of field 
strengths realized over the visible part of the surface of the 
white dwarf is given by the top grey bar. The left-hand 



panel shows the B-tp diagram, a greyscale plot of the fre- 
quency distribution of the magnetic field strength B and 
the direction cosine cosip. The fractional contribution of 
each single database spectrum to the integrated spectrum 
is represented by the greyscale value of the corresponding 
pixel in the plot. This presentation includes the effects of 
pixel area, foreshortening, and limb darkening. The sum 
of all filling factors would be unity if limb darkening were 
neglected, but falls below unity with limb darkening in- 
cluded. The B-t/j diagrams are equivalent to the ZEBRA 



plots of Donati et al. (1994), except for the effect of limb 
darkening which was not included by these authors. The 
diagrams illustrate the change in the weighting of the two 
main database parameters, B and ip, as the star rotates. 
The right-hand panel shows the resulting integrated flux 
and circular polarization spectra at the four rotational 
phases. 

Model (Al), centred dipole viewed at i = 60°: The po- 
lar field strength is £? d ol = 40 MG and the axis points to- 
wards (9 d , $ d ) = (60°, 340°). This oblique rotator model 
stands for a simple low-field geometry. The hidden fraction 
is only 7% of the white dwarf surface. The flux spectra in 
Fig. [|c are quite similar at <f> = 0.25, 0.5, and 0.75, but the 
circular polarization spectra are not. The B— tp diagram 
looks different at <fi — and so does the flux spectrum. 
These differences suggest that full phase coverage is es- 
sential for a successful recovery of the field geometry. 

Model (A2), centred dipole viewed at i = 30°: The hid- 
den fraction of the surface is now 25%. Otherwise, the 
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Fig. 5. Top: Model (A2), centred dipole viewed at i = 30°: a) B-ip diagram, b) field strength and longitudinal 
component, c) flux and polarization spectra. Bottom: Model (B), pure quadrupole viewed at i = 60°: a) B— ip diagram, 
b) field strength and longitudinal component, c) flux and polarization spectra. See Fig. for further explanation. 



properties of the model (Fig. [5], top) are similar to (Al). Model (B), centred quadrupole viewed at i = 60°: The 

For the centred dipoles of models (Al) and (A2), the cir- polar field strength is -Bp ol = 40 MG and the axis points 

cular polarization vanishes at a phase <pQ, where the dipole towards (0 q , <& q ) = (60°, 340°). Figure ||c (bottom) shows 

axis is oriented perpendicular to the line of sight. that there is little rotational variation. The flux spectra 
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Fig. 6. Top: Model (C), aligned dipole and quadrupole viewed at i = 60°: a) B— tp diagram, b) field strength and 
longitudinal component, c) flux and polarization spectra. Bottom: Model (D), shifted high-held dipole viewed at 
i = 60°: a) B—ip diagram, b) field strength and longitudinal component, c) flux and polarization spectra. See Fig. |4| 
for further explanation. 



and the polarization vary little for cj> = 0.25, 0.5, and 0.75, 
but differ at = 0. 

Model ( C), aligned centred dipole and quadrupole 
viewed ati = 60° : A quadrupole of SjJ , = 20 MG is added 



to a dipole of B d ol = 40 MG with (9 d , $ d ) = (9 q , $ q ) = 
(60°, 340°). The asymmetry introduced into the field ge- 
ometry causes larger rotational variations in flux and po- 
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Fig. 7. Top: Model (E), shifted dipole viewed at i = 60°: a) B-ip diagram, b) field strength and longitudinal com- 
ponent, c) flux and polarization spectra. Bottom: Model (F), non-aligned dipole-quadrupole combination viewed at 
i = 60°: a) B-ip diagram, b) field strength and longitudinal component, c) flux and polarization spectra. See Fig. |^ 
for further explanation. 



larization than for the pure dipole or quadrupole (Fig. |gc, 
top). 

Model (D), shifted high- field dipole viewed at i = 60°: 
The polar field strength is B* ol = 110 MG, offset in all 



three coordinates by (x' s, y' Q g, z'off) = (0.05, —0.10, 
0.15). The shift along the dipole axis increases the maxi- 
mum field, and the sideways shift decreases the minimum 
field to the effect that B ranges from 39 to 192 MG. That 
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is, B varies by a factor of five compared to a factor of 
two for the centred dipole. The high field causes the flux 
spectra to show substantially less structure than in the 
previous models, suggesting that a higher signal-to-noise 
ratio is needed for reconstruction (Fig. |6|c, bottom). There 
is substantial variation in the circular polarization over the 
rotational period, however, which helps in the reconstruc- 
tion. 

Model (E), shifted dipole viewed at i = 60°: This is 
an extremely off-centred dipole with J3 d ol = 58 MG and 
(x'ofi, 2/ off; z 'otf) = (0.15, —0.10, 0.30) which displays a 
variation of B over the surface by nearly a factor of 12. 
For one half of the rotational period, the high-field pole is 
in view, over the other half the field distribution is concen- 
trated at low field strengths. Effectively, this represents a 
star with a low field of around 20 MG over most of the star 
and a spot in which the field rises to 198 MG. The circu- 
lar polarization displays pronounced rotational structure 
(Fig. 0c, top). 

Model (F), non-aligned dipole- quadrupole combination 
viewed at i = 60°: This is the most complex field model 
featuring the superposition of a non-aligned dipole and 
quadrupole with equal polar field strengths of 40 MG. The 
polar directions, (9 d , $ d ) = (60°, 340°) and (9 q , $ q ) 
= (30°, 250°), are separated by 64°. The field geometry 
features two high- field spots, an upper positive and a lower 
negative one, which are dominated by the quadrupole and 
the dipole, respectively, and are separated by ~90° (Fig. [j], 
bottom). 

An overview of the models (A)-(F) is given in Table § 
To illustrate the effects of noise in the spectra used in 
the reconstructions, we show in Fig. || the flux and polar- 
ization spectra of Model (Al) for cf> = at noise levels of 
S/N = 100 and S/N = 20. 

5.2. Optimization algorithm 

Our spectral synthesis method is sufficiently fast to allow 
the use of hierarchical search strategies in the parameter 
space. We utilize the optimization routine evoC (Trint & 
Utccht 1994[) that impleme nts an evolutionary strategy al- 
gorithm (Rechenberg 1994), and has proven useful already 
in other astrophysical contexts (Gansicke & Beuermann 
19961; |Gansicke et al. 1998j; [Kube et al. 2000[). The tasTTs 



to find a set a = (ai, . . . , ajv/) of M free parameters that 
minimizes the classic penalty function 



P 3 



on 






-0.! 
0.1 

°l 
0.1 

4000 



S/N=10Q 



H 1 1 H 1 1 1 1 1 1 1 1 1 1 1 1 1— 



, 1 1 1 , 1 1 1 1 , 1 1 1 , 

5000 6000 7000 

Wavelength / k 



1 



JY 



M - M 



s 



LL 



doYv 



Fig. 8. Examples of simulated spectra used as input to the 
reconstruction procedure. The cf> = spectrum of Model 
(Al) is shown with a signal-to- noise ratio S/N = 20 and 
100 (shifted upwards by one flux unit). Top: Flux spectra, 
bottom: circular polarization spectra. 



For each field model and each reconstruction with a 
certain set of free parameters, the evoC optimization pro- 
cess is run repeatedly, typically 6-20 times, starting each 
run with different, randomly chosen parameter values. 
Not all runs end up in the global minimum. A misguided 
run may be caught in a local minimum corresponding to 
an incorrect field configuration, which nevertheless has a 
Zeeman spectrum similar to the input one. We define a 
success rate of the optimization as the fractional number 
of runs which reach a best-fit x? e d < 2 -° for S/N = 100, 
and Xred < 1-1 f° r S/N = 20 (corresponding virtually al- 
ways to the global minimum). As a last finish, we employ 
a downhill simplex algorithm ( Nclder & Mead 1965|; |Presl 



given the input data pixels fj, the model data pixels s 



et al. 1992[) on the run with the best xjr cd , which sometimes 



improves on the evoC solution. 

In order to illustrate the problem associated with lo- 
and the standard deviations a). Good fits require xL^ 1 cal minima in the X 2 -landscape, we present in Fig. | a 



j ' 



We have applied Gaussian noise to the input spectra to 
yield signal-to-noise ratios, corresponding to relative stan- 
dard deviations in F\ and absolute standard deviations in 
V/I, of 0.01 and 0.05, respectively. For fits to fluxes only, 
j runs up to N = 1444, and, for fits to both flux and po- 
larization, up to N — 2888. 



simple example of different field geometries which yield 
similar Zeeman spectra. The input geometry is the sum 
of a (non- rotating) dipole with B^ ol = 40 MG and an 
aligned quadrupole with B^ ol = 20 MG, viewed at i = 60° 
(with 9 d = 9 q = uj = 0°, Fig. |). The right-hand 
panel of Fig. U shows the corresponding flux and polar- 
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Fig. 9. Left panel: Contour plot of the x 2 -landscape for the spectral flux in the £?p ol ,0 q -plane. The input configuration 
is given by B^ ol = 40MG, £? q ol = 20 MG, and 8 d = 9 q = 0°. Darker shading corresponds to smaller values of x 2 c d- 
Centre panel: B-ip diagrams of the field configurations corresponding to the local (top) and the global (bottom) 
minimum. Right panel: Flux and circular polarization spectra corresponding to the the global minimum (lower curves) 
and the local minimum in the upper left of the left panel (upper curves, shifted upwards by 0.5 units in flux and 0.1 
units in polarization). 



ization spectra (lower curves). We add Gaussian noise 
of S/N = 100 and compute flux and flux+polarization 
spectra covering a range of quadrupole parameters, with 
the dipole parameters and the inclination kept fixed. The 
quadrupole is allowed to vary in strength and orientation 
with <d q free at <E> q = 90°. For this choice of parameters, 
the quadrupole is perpendicular to the dipole and to the 
line of sight for 6 q = 90°. The left-hand panel of Fig. | 
shows a contour plot of the x 2 -landscape for the spec- 
tral flux in the -Bp ol ,9 q -plane. Besides the global mini- 
mum at the parameter values of the input configuration 
(_Bp j = 20 MG, 8 q = 0°), a second pronounced minimum 
appears at B^ ol = -17MG and 9 q = 90°, with the minus 
sign indicating a reversed polarity of the quadrupole. The 
flux and circular polarization spectra for this minimum are 
also shown in the right-hand panel of Fig. ^. At moder- 
ate noise levels, the flux and polarization spectra of these 
two diverse field geometries become indistinguishable and 
it is not surprising that the local minimum (in the up- 
per left corner of the contour plot) persists if flux and 
circular polarization are considered together. The shallow 
local minimum at 5 q ol ~ 33 MG, 9 q ~ 80°, on the other 
hand, disappears if V/I is included in the computation 
of Xred- I n tne cen t re panels of Fig. [|, the B-vp diagrams 
for both configurations are shown. Both distributions are 
sufficiently similar if projected either on the i?-axis or on 
the cos ^-axis to explain why the spectra are similar, but 
not identical. 

Finally, we note that fitting the remaining parameters 
of the field model (like B^ ol ) instead of keeping them fixed 
would cause the local minima to become even more pro- 
nounced. Increased noise also deepens the local minima 
relative to the global one. A local x 2 -minimum is respon- 



sible for an incorrect, although not entirely dissimilar, re- 
construction of Model F discussed below. 

5.3. Reconstruction fits 

5.3.1. General characteristics of the solutions 

Depending on the complexity of the input field, we con- 
sider reconstructions which differ in the numbers of free 
parameters, ranging from the full set of ten down to seven 
(with the quadrupole component or the offset neglected). 
Some redundancy is allowed because a dipole offset along 
its axis can also be modelled, to first order, by an aligned 
centred dipole-quadrupole combination. With data of suf- 
ficient S/N, the reconstruction procedure can recognize 
such subtle differences. 

As a general feature, the reconstructed global field is 
of relevance only for that part of the stellar surface which 
is visible during the observation (or covered by the syn- 
thetic input in this paper). This underlines the importance 
of phase-resolved observations which allow the determina- 
tion of the inclination and, thereby, to estimate the oc- 
culted fraction of the surface. 

Spectrophotometry of high S/N is obtained more eas- 
ily than spectropolarimetry of the same quality. The ob- 
server has to decide, therefore, whether a given amount of 
observing time is better spent on high-quality spectropho- 
tometry or on circular spectropolarimetry of lower qual- 
ity. In order to address such questions, we reconstructed 
all field geometries, using the spectral flux and the po- 
larization, and using the spectral flux only. The flux-only 
reconstructions are successful in several cases, but the de- 
viations from the input geometry tend to be larger, and 
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an increased number of non-convergent runs suggest a less 
well-behaved x 2 -landscape. We find that the circular po- 
larization is not needed in simple cases, while its inclusion 
is extremely useful for the reconstruction of more complex 
fields. 

5.3.2. Results for individual field geometries 

In this Section, we present the results for the reconstruc- 
tions of the input models (A) to (F), using the spectra 
at four rotational phases. All results are listed in Table 1. 
The column denoted 'flag' indicates whether the fit is to 
flux and polarization (fp) or to the flux only (f). The 
last column illustrates the convergence properties in the 
multidimensional parameter space, referred to as success 
rate above (number of successful runs vs. total number of 
runs). 

Model (Al), centred dipole viewed at i = 60°: The re- 
sults in lines 1-4 assume B^ ol = 0, those in lines 5-8 zero 
offset. All reconstructions are successful and reproduce 
the dipole field strength within 0.1 MG and the magnetic 
axis and the inclination with rms deviations of 5°and 8°, 
respectively. Not surprisingly, the accuracy of the recon- 
struction benefits from a high S/N, but is acceptable even 
for flux-only fits and a low S/N ratio. Note that errors in 
$ d are irrelevant as long as d matches closely. The same 
holds for <& q and 9 q as long as B^ ol is close to zero. If 
all parameters are included in the fit (lines 9 and 10), 
a quadrupole component usually appears which is largely 
compensated for by a shift in the dipole (plus quadrupole) 
to the effect that the net field is dipole- like again. The low- 
noise flux-and-polarization fit of line 9 is quite acceptable, 
while the high- noise flux-only fit of line 10 produces larger 
misfits in i and in the field geometry. Even the latter pro- 
vides an acceptable reconstruction over the visible part of 
the surface, but deviates strongly from the input in the 
permanently occulted part. This result is due to the in- 
clusion of a higher multipole component than present in 
the input. Figure [l^ demonstrates this result. 

Model (A2), centred dipole viewed at i = 30°: The re- 
constructions are of a quality similar to that of model 
(Al). The permanently occulted fraction of the stellar sur- 
face, for which the reconstruction remains undefined, is 
now 25%. 

Models (B, C), pure quadrupole and aligned dipole- 
quadrupole viewed at i = 60°: The four reconstructions 
each use the same set of parameters as lines 5-8 of model 
(Al). The lack of a dipole component in (B) and the rel- 
ative strength of the dipole and quadrupole components 
in (C) are recognized even in the flux-only and high-noise 
fits. 

Model (D), shifted high- field dipole viewed at i = 60°: 
The dipole is aligned with the rotation axis. Hence, the 
rotational modulation is caused by the off-centre shift. 
In spite of the weaker Zeeman structures (Fig. ^c, bot- 
tom), the configuration is recovered correctly if the search 
assumes a shifted dipole (lines 1-4). The absence of a 




Fig. 10. Aitoff maps of the magnetic field strength B (in 
MG) for Model (Al). The input (top) and the high-noise 
flux-only reconstruction of Table [|, Model (Al), line 10 
(bottom) are shown. The axis of the maps represents the 
rotation axis. The centre of the maps marks the —xo di- 
rection (see Fig. |^b). The region within 30° from the lower 
pole is permanently hidden (dashed white line). 



quadrupole component is recognized in the low-noise flux- 
and-polarization fits (line 5), but less so in the high- noise 
flux-only fit (line 6). Interestingly, the offset is recognized 
correctly in both cases. 

Model (E), shifted dipole viewed at i = 60°: This is the 
model which features one high-field spot. If modelled as 
a shifted dipole, the parameters are recovered with high 
accuracy, even when the S/N is low and the polarization 
is disregarded (lines 1-4) . Allowing for a quadrupole com- 
ponent leads to the usual compensatory effects (lines 5 
and 6). The reconstruction is acceptable over the visible 
surface, and deviates only slightly from the input field in 
the occulted part. 

Model (F), non-aligned dipole- quadrupole combination 
viewed at i = 60° : All fits reproduce the general structure 
of the field with its two spots, separated by ~90°, but only 
the fit using flux and polarization spectra at S/N = 100 
(line 1) correctly finds the axes of both, the dipole and the 
quadrupole components. Convergence problems, which 
arise when the polarization is disregarded (line 3), may 
be due to a more corrugated x 2_ landscape compared with 
the fits including the polarization. Figure ^ provides an 
overview of the results for model (F). The input is depicted 
in Fig. pT| a and the line-1 reconstruction in Fig. pT|b. The 
high- noise flux-only fits of lines 4 and 6 (Figs. |TT| c and|Tl"|d) 
deviate in the strengths and orientations of the dipoles and 
quadrupoles. In Figs. |ll|a and |TT|b, the high-field spot on 
the northern (southern) hemisphere is dominated by the 
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Fig. 11. Aitoff maps of the field strength B (in MG) showing the results of reconstruction fits for the Model (F) 
field distribution (non-aligned dipole-quadrupole combination). The centre of the maps marks the —xq direction. The 
axis of the quadrupole and the field direction at its poles is indicated by the + and — symbols, the axis of the 
dipole correspondingly by the circled + and — symbols (see text), a) Input field, b) low-noise flux-and-polarization 
reconstruction of Table |l|, Model (F), line 1, c) high- noise flux-only reconstruction of line 4, and d) high- noise flux-only 
reconstruction of line 6. 



quadrupolar (dipolar) contribution. In Fig. |li"|c, the po- 
larity of the quadrupole is reversed. Finally, in Fig. |ll|d, 
dipole and quadrupole are nearly aligned, but the spots 
are shifted due to the finite values of x' ff and y' oS . 
Nevertheless, these fits are not altogether wrong if only 
the two-spot structure is considered. They suggest that 
a low Xred does, in fact, indicate a representation which 
bears some similarity to the input field structure, even 
if the choice of multipoles differs from that of the input. 
Seemingly, better results cannot be expected given the 
high noise of these two fits. 



In summary, the code is able to reconstruct the mag- 
netic field geometries of the type discussed here from 
phase-resolved flux and polarization spectra of high S/N. 
Our experience is that the inclusion of more than four 
phases does not improve the fits substantially, which is 
understandable given the overlap in surface coverage. 
Naturally, the reconstruction becomes less perfect when 
the polarization information is excluded and the noise is 
increased. Leaving off the polarization also seems to cre- 
ate convergence problems. An important aspect is that 
phase-dependent spectra allow the inclination of the rota- 
tion axis to be determined along with the field geometry. 
The accuracy achieved ranges from a few degrees to about 
20° depending on the quality of the spectra. 



5.3.3. Fits to a single spectrum 

If only a single set of flux and polarization spectra or a 
single flux spectrum is available, information on the field 
distribution is reduced to the visible hemisphere. The lo- 
cation of the rotation axis remains unknown, and only the 
angle between the magnetic axis and the line of sight is 
constrained by the fit. We have performed similar tests as 
above to single sets of spectra and find that simple field 
geometries can still be recovered. 

6. Discussion 

We have presented a formalized approach to the inter- 
pretation of phase-resolved flux and circular polariza- 
tion spectra of rotating magnetic white dwarfs (MWDs). 
Tomographically locating positions with a certain field 
strength B on the surface of the star is hampered by the 
fact that only the self-eclipse of such a region manifests 
itself in flux spectra, while the positional information con- 
tributed by the rotation is obliterated by the Stark broad- 
ening. 

6.1. Present approach 

In our approach, we determine the parameters of a global 
field model directly by a least-squares fit to the spec- 
tral data. We caution that it is not a priori clear to 
what extent the global field can be constrained by such 
an approach, because the spectral information represents 
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an average over the visible hemisphere at each phase. 
Our results demonstrate, however, that the phase-resolved 
Zeeman spectra contain enough information to allow the 
reconstruction of the field geometries considered by us. 
These involve combinations of dipoles and quadrupoles 
which are allowed to have different axes and to be shifted 
off-centre. The model contains up to ten free parame- 
ters and is sufficiently general to allow for rather com- 
plex surface field geometries featuring, e.g., a dominant 
single high-field spot, two spots separated by much less 
than 180°, or even a bipolar spot on an otherwise low- 
field star (not included in the models presented here). An 
advantage of our approach is that these fields automati- 
cally fulfil the requirement of being produced by sources 
inside the star. A disadvantage is the limitation in the 
number of free parameters. 

In addition to the cases presented here, we have also 
attempted to reconstruct octupolar fields and were suc- 
cessful for aligned dipole-quadrupole-octupole combina- 
tions. However, if all multipole components with I = 3 
and m = 0, . . . , 3 are included (15 parameters for the ex- 
pansion, two angles describing the direction of the refer- 
ence axis, and the inclination), the evoC minimization al- 
gorithm encounters convergence problems, caused by too 
large a number of free parameters. 

Another important aspect is the level of the signal-to- 
noise ratio S/N required for a successful reconstruction of 
the field. The model atmospheres of hydrogen-rich MWDs 
are characterized by the rather broad and strong Zeeman- 
shifted Balmer lines which allow a field reconstruction 
already for S/N = 20-100: S/N = 20 is the lower limit, 
while there is little improvement for S/N > 100. For com- 
parison, Zccman-Doppler imaging of main sequence stars 
operates on much fainter meta l lines and needs a much 
higher S/N ( Brown ct al. 1991 ). However, while Zeeman- 
Doppler imaging is performed over individual lines, the 
high field strengths of the MWDs require a fit over the 
whole visible wavelength range. 

In the analysis of observed Zeeman spectra, one may 
encounter some problems which are absent in the present 
reconstruction of synthetic spectra. While the variation 
of the statistical noise amplitude with wavelength can be 
accounted for in the x 2 -statistic (Eq. |l3|), systematic un- 
certainties between the observed and calculated spectra 
cannot: (i) errors in the theoretical database spectra; and 
(ii) errors in the flux calibration of the observed spectra. 
Errors of type (i) may prevent a satisfactory convergence 
of the fits and/or lead to incorrect values of the parameters 
describing the field. Our experience is that such errors are 
of minor importance, given the present state of the the- 
ory of radiative transfer in magnetic stellar atmospheres. 
Errors of type (ii) may affect the ability to recognize high 
field spots on stars with a predominantly moderate field. 
For example, in model (F) at phase <fi — 0.25, and simi- 
larly in other models, the Ha a~ component consists of 
a shallow depression extending from 5000 to 6000 A. An 
error in the flux calibration which happens to weaken or 
strengthen such a depression can lead to serious errors in 



the derived field distribution. A careful flux calibration is, 
therefore, of utmost importance. 

6.2. Different optimization strategies 

Should MWDs turn out to have field geometries which 
are more complex and require more free parameters than 
adopted by us, we may have to consider alternative op- 
timization techniques. E.g., a genetic algorithm may be 
more robust than the evoC code and allow for a somewhat 
larger number of parameters. A full harmonic expansion 
with I 3> 2 may become tractable if a regularization oper- 
ator like MEM drives the solution towards low-order fields 
while permitting higher-order components to be used as 
necessary to fit the data. 



6.3. The ZEBRA approach 

The problem of retrieving the field structure of rotating 
MWDs has previously been studied by Donati et al. (1994) 
in what they called the ZEBRA approach. They used a 
maximum-entropy method to deduce the most likely two- 
dimensional frequency distribution f(B t ,B\) of the trans- 
verse and longitudinal field components with respect to 
the line of sight separately for each rotational phase. The 
method has the obvious advantage that no a priori as- 
sumption is made about the global field structure. On the 
other hand, the interrelation between the overlapping field 
distributions at different rotational phases is not utilized 
and there is no prescription for the interpretation of such 
an interrelation in terms of a global field. Indeed, there 
is no guarantee of a physically meaningful reconstruction 
(e.g. sources only within the star leading to a curl-free 
field outside the star). Thus, the detailed structure and 
the physical characteristics of the global field rem ain un- 
defined in t he ZEBRA method in the present form ( Donati 
|et al. 1994| ). 

If the underlying global field structure is sufficiently 
simple, it may be derived in a second step added to the 
ZEBRA method. In a first step, the best-fit ZEBRA dia- 
grams (or B-?p diagrams similar to ours) are determined 
using a MEM-type regularization scheme as suggested by 
Donati et al. (1994). In a second step, a parametrized 



global magnetic field model is then fitted to the phase- 
resolved ZEBRA diagrams. Since the second step would 
not involve the computation of spectra from the database, 
which is by far the most time-critical process in the present 
method, this two-step approach is probably advantageous 
with respect to CPU time. Without detailed tests, how- 
ever, it is not clear whether this approach would be supe- 
rior to directly fitting the Zeeman spectra. 

6.4. A future approach 

One may endeavour to relax the restrictions on the 
global field structure by parametrizing the surface field 
as (Bi, . . . , -B_/v) for a star with N surface elements and 
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to impose a regularization scheme, e.g. MEM, to ensure 
the smoothness of the solution. The feasibility of such an 
approach, its convergence properties, and the interpreta- 
tion of the derived field model would still have to be stud- 
ied, however, as well as the demands on computation time 
given the formidable number of parameters. 

7. Conclusion 

We have described a method to reconstruct the field struc- 
ture of magnetic white dwarfs which provides an internally 
consistent fit to spcctropolarimctric data taken at different 
rotational phases in terms of a parametrized field model. 
We presently use dipoles and quadrupoles which are al- 
lowed to be shifted off-centre to increase the versatility of 
the model. An application to real data will be described 
elsewhere. 

We do not know whether MWDs have the regular fields 
adopted here or possibly field structures as complex as 
spotted main sequence stars. Fortunately, there are several 
single white dwarfs with known rotational periods, and 
about 60 rotating MWDs in cataclysmic binaries, some of 
which are known to have fields which deviate from sim- 
ple centred dipoles. The study of such systems using the 
present and similar techniques promises to increase our 
knowledge of the end-product of magnetic stellar evolu- 
tion. 
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Table 1. Reconstructed magnetic parameters for the configurations (A)-(F). Each model is introduced by a boldface 
line which lists the input parameters. The subsequent numbered lines represent the individual reconstructions. In the 
'flag' column, 'fp' denotes simultaneous fits to flux and polarization spectra, 'f fits to flux spectra only. The last 
column indicates the success rate of the convergence of the multidimensional parameter search. Note that for (B) and 
(C), 6^ = 6 d and = $ d . 
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